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Abstract - We address the problems in applying cycle expansions to bound chaotic 
systems, caused by e.g. intermittency and incompleteness of the symbolic dynam- 
ics. We discuss zeta functions associated with weighted evolution operators and in 
particular a one-parameter family of weights relevant for the calculation of classical 
resonance spectra, semiclassical spectra and topological entropy. For bound intermit- 
tent system we discuss an approximation of the zeta function in terms of probabilities 
rather than cycle instabilities. This approximation provides a generalization of the 
fundamental part of a cycle expansion for a finite subshift symbolic dynamics. This 
approach is particularly suitable for determining asymptotic properties of periodic 
orbits which are essential for scrutinizing the semiclassical limit of Gutzwiller's semi- 
classical trace formula. The Sinai billiard is used as model system. In particular 
we develope a crude approximation of the semiclassical zeta function which turns 
out to possess non analytical features. We also discuss the contribution to the semi- 
classical level density from the neutral orbits. Finally we discuss implications of our 
findings for the spectral form factor and compute the asymptotic behaviour of the 
spectral rigidity. The result is found to be consistent with exact quantum mechanical 
calculations. 



1 Introduction - Dilemmas in Quantum Chaology 

There is no definition of Quantum Chaology or Quantum Chaos which is generally agreed upon. 
One may focuse on the correspondance principle and state the question How is chaos revealed 
in a quantum system when the Planck's constant h —^ 0, if its classical counterpart is chaotic. 
The answer turns out to as complicated as the question is interesting. It seems as if one cannot 
understand properly the transition from quantum to classical mechanics for a system with few 
degrees of freedom without considering its coupling to the environment. This coupling may induce 
decoherence, enabling quasi classical behaviour [Q . 

One may also (and we will) focuse on the more modest question: What is the behaviour 
of a quantum system if its classical counterpart is chaotic? We thus avoid any discussion of 
decoherence phenomena. Restricting our attention to bounded autonomous systems we may 

^Proceedings of the Los Alamos Center for Nonlinear Science Quantum Complexity in Mesoscopic 
Systems, May 1994, to appear in Physica D. 
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choose e.g. to study quantum spectra. There has been a lot of numerical and experimental 
studies of the statistical properties of quantum spectra of classically chaotic systems. There 
are some widely known conjectures inspired by random matrix theories [||. Although some of 
the numeric evidence is rather striking there is very little theoretical understanding. Such an 
attempt naturally starts from the semiclassical trace formula. This formula relates the density of 
eigen-energies to dynamical invariants of the periodic orbits of a system Q : 

Here g{E) — Tr G{q,q'; E) is the trace of the semiclassical Green's function, and the density of 
eigenstates is given by d{E) = — ^ limg^o iJxig{E + ie). The index p labels the primitive periodic 
orbits, Sp is the action integral along the orbit, Tp its period, Mp is the linearized Poincare map 
around the orbit (monodromy matrix), jj,p is the Maslov index, and go{E) provides the mean 
level distribution. We will exclusively consider billiards and put h = m = 1. The semiclassical 
limit ?i — > is thus replaced by -B — s- cx). The action is given by Sp — Ip ■ k{E) where k is the 
momentum k — V^E, and Ip is the length of the prime orbit. 

The trace formula is derived through a stationary phase approximation, except for goiE) 
which is calculated from the short time behaviour of the semiclassical propagator. In principle 
goiE) is given by an asymptotic expansion whose leading term is the Weyl term. The division 
into a mean contribution and an oscillating may turn out to be more involved than anticipated. 
We will indeed encounter such problems in this paper. 

It is often more convenient to study the Gutzwiller-Voros zeta function Q which for systems 
with two degrees of freedom reads 

where Ap is the expanding eigevalue of Mp. Zqy is related to the level density by 

g{E)-g^{E)^^\uZ{E) , (3) 

so the zeros of the zeta function are the poles of the of eq. ([|), i.e. the semiclassical eigenvalues. 

It is not known whether the trace formula manages to give the right density of states on, 
or close to the real energy axis. Whatever method one uses, only a moderate number of zeros 
can be computed from a reasonable number of periodic orbits. The exponential proliferation 
of periodic orbits makes every step towards the semiclassical limit very expensive. It therefore 
seems impossible in practice to study the semiclassical limit of the semiclassical trace formula by 
explicit computation of the periodic orbits. 

This motivates a study of asymptotic measures of the set of periodic orbits for bound systems, 
which is the main problem to be addressed in this talk. 

We will make extensive use of the concept of zeta functions. The most succesful use of zeta 
functions has been achieved for so called Axiom-A systems (a very strong definition of chaos 
suitable for making mathematical theorems). Then there are theorems about entireness of the 
certain classes of zeta functions. However, bounded systems are generally not Axiom-A, they are 
generally intermittent and lack a simple symbolic dynamics, so practical applications of the theory 
of zeta functions are not very successful. Since our main motivation is to understand fluctations 
in the spectra of bound systems we are presented with another dilemma. It will however turn out 
that we can make use of this inconvenient intermittent property to formulate approximate zeta 
functions. 



2 A menagerie of zeta functions 

This exposition of zeta functions and the thermodynamic formalism is rather brief. The reader 
may found more details in refs ^, 0, ^, ^ and references therein. 
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2.1 Evolution operators, trace formulas and zeta functions 

Zeta functions are introduced by considering the evolution operator . It describes the evolution 
of a phase space density <I>(x) 

Cl^{x) = J wix,t)S{x ~ fiy)My)dy . (4) 

The phase space point x is taken by the flow to /*(a;) during time t. w{x,t) is a weight 
associated with a trajectory starting at x and evolved during time t. It is multiplicative along 
the flow, that is w{x,ti + ^2) = ti)w(/*i (x), ^2)- This ensures that the eigenvalues are on 
the form A^. 

As we are only studying billiards it is convenient to use the length traversed I as time variable. 
We now compute the trace of the evolution operator, that is, the sum of its eigenvalues 
trCl = J2a e*''"'- The trace may be written as a sum over the periodic orbits in the system 

trC'^ ^ J w{x,t)S{x - f{x))dx = 
^ djl-nlp) 

p n—l ' ^ P ^ ' 

The trace may also be written as the Fourier transform of the logarithmic derivative of a zeta 
function 



The zeta function reads 



p m— 



^.(fc)=nn ' (7) 



The ka 's introduced above should be the zeros of the zeta function. For the system we are going 
to consider things will be complicated and the operator will not have a discrete spectrum. 

By using different weights w one can probe different properties of the flow. We will consider 
the one-parameter family of weights uo = \K{x^t)Y . K{x,t) is the expanding eigenvalue of the 
Jacobian transverse to the flow. It is only approximately multiplicative along the flow but it 
is possible to modify it slightly so as to become exactly multiplicative [|l^. However, this is a 
subtlety compared to the approximations we are going to apply. 

The leading zero will generally be on the imaginary axis k = —ih{T), where /i(t) are gener- 
alized entropies, in the thermodynamic formalism they are usually named topological pressures, 
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Putting T = we obtain the classical zeta function whose zeros yields the so called resonance 
spectrum or spectrum of correlation exponents (again things are complicated if the spectrum is 
not discrete) 

The leading zero provides the escape rate. For a bound system the escape rate equals zero: 
h{0) = 0. 

The topological zeta function is obtained by considering the case t — 1. The leading zero now 
gives the topological entropy h{l). The asymptotic behaviour of the trace is 

00 

trClj, ^J2lpJ25{l- nip) ^ e'^d)* (9) 

p n—l 



3 



so that the number of cycles with periods less than I is ~ e'''^-''//i(l)/. To obtain this exponential 
increase, the contour of the Fourier transform must extend below all zeros where the Euler 
product (0) converges. 

Let us now study the case r — 1/2. We then get 

This zeta function is called the quantum Fredholm determinant Q and is equivalent to the 
Gutzwiller-Voros zeta function in the semiclassical limit through the identification k{E) — —k. 
Its leading zero is closer to the origin than the topological one < h{l/2) < h{l), cf. sections 
2.2 and 3.3. This zero provides a pole in the trace formula, and the periodic orbit sum (|^) as 
well as the trace formula ([l]) diverge in the half plane above it. This is usually referred to as a 
an entropy barrier. The Euler-product representation of the zeta function ^ also diverges above 
this leading zero (because the zero is not a zero of an individual factor and an infinite product 
is not allowed to converge towards zero). We have left out the Maslov indices but it is possible 
to account for them in the weight as well. The introduction of phase indices will shift the the 
leading zero towards the origin as we soon will see. 



2.2 Cycle expansions of zeta functions 

The zeta functions have been formulated as Euler products. These products diverge wherever 
they are interesting (i.e. at the nontrivial zeros). In computations, one generally expands them 
into infinite series, power series in the case of maps and Dirichlet series in case of flows. Power 
series converge up to its first singularity and are thus better for computations. Dirichlet series 
are more complicated, they may have a strip of conditional convergence and do not neccesarily 
have singularities on its border(s) of convergence. 

Zeta functions are entire if the system has a Poincare map which is analytic, hyperbolic and 
have a finite subshift symbolic dynamics ]l^ . In that case the expansion may be divided into 
a fundamental part, yielding the gross structure of the spectrum, and curvature corrections [^j. 
The standard example is the expansion of the dynamical zeta function (i.e. the m = factor of 
the zeta functions above) for a system with binary symbolic dynamics: 

np(i - 1,) = 

l-ti-to 
^[tio ^ tito] 
^[^100 ^ ^lo^o] - 
— [^1100 ~ titioo 



\f tf ] (1^) 

- iiio^o + titioto] ~ [iiiio — titiio] — [iiooo — ^loo^o] 



where tp — exp(— i^pfc)/|A)3|^^'^. The expansion is organized according to increasing symbol 
length in an obvious way, in the case of a map (the Zp's being integers) this means according to 
increasing powers of exp{—ik). The expansion is dominated by the fundamental contribution, 
1 — ti ~ to. For the lowest powers there is a nice interpretation why the rest of the expansion is 
small 1^] . The tail of the expansion is organized in such a way that each square bracket contains 
some long orbit(s) (length n) minus its approximant(s) in terms of shorter ones. The sizes of 
these curvature corrections fall of exponentially with n. This leads to a pole lying beyond the 
leading zero. Similar expansion for the full zeta function has infinite radius of convergence (under 
the above conditions) such that nonleading zeros may be extracted. 

Example. Let us consider a hypothetical billiard having the Baker's map as Poincare map. 
Then |Ap| = 2"p. Let further the periods be discretized Tp — UpTg. Only the fundamental part 
will remain: 

gTofc 
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with leading zero = — irlog2/To. We see that h{0) = as it should be for a bound system and 
the topological entropy is h{l) = log2/ro- This is of course a very unrealistic example and the 
semiclassical spectrum, t — 1/2, does not look sensible at all. 

If the symbolic dynamics is a more general, finite subshift (meaning that there is a finite set 
of pruning rules), this scheme may still be worked out. There is a wealth of illustrative examples 
in ref. However, if the symbolic dynamics is not a finite subshift there will be no similar 
expansion. The scheme outlined here is therefore not applicable for a generic bound system. In 
this case there is not even a well defined fundamental part. We will suggest a solution of this 
dilemma in the next section. 

2.3 Zeta functions in the BER approximation 

In ref. an approximate expression for the zeta function is given for intermittent, ergodic 
Hamiltonian systems. The idea is based on a paper by Baladi, Eckmann and Ruelle so we 
refer to it as the BER approximation. 

In an intermittent system laminar intervals are interupted by chaotic outbursts. Let 
be the time ellapsed between two succesive entries into the laminar phase. The index i labels 
the i'th interval. Provided the chaotic phase is chaotic enough, the lengths of the intervals 
are presumed uncorrelated, and A may be considered as a stochastic variable with probability 
distribution p(A). The zeta functions may then be expressed in terms of the Fourier transform 



We use the symbol to mark all quantities derived in the BER approximation. 

3 The Sinai billiard - Classical considerations 

The Sinai billiard |jlj| consists of a unit square with a scattering disk, having radius R; < 2R < 1, 
centered on its midpoint. This billiard has a fairly simple geometry, but exhibits many features 
typical for bound chaotic system; it is intermittent and lacks a simple symbolic dynamics. 

The trajectory of a particle in the Sinai billiard consists of laminar intervals, (bouncing 
between the straight sections) interrupted by scatterings off the central disk. The Sinai billiard 
seems therefore suited for the BER approximation. The variable A introduced above is simply 
the length of the trajectory segment between two disk bounces. 

The material of this section, except for sect. 3.3, is developed in ref in more detail. 

3.1 Formulation of p(A) and a generalization 



We use the disk to define the surface of section. The canonical variables are {R<j), V2E sin a) with 

the angles < (j> < 2tt and — ^ < a < ^ defined in fig. la. The map (0, sina) t-^ (0, sina) is 
area preserving and has uniform invariant density. 

We will now use this uniformity to formulate an expression for p(A). To that end we use 
a natural partition of the surface of section defined by the topology of the problem. Suppose a 
laminar segments hits the vertical walls Ux times, and the horizontal walls riy times. Equivalently 
one can unfold the billiard into a regular lattice of disks. These disks are assigned numbers 
q — {nx,ny), as in fig lb. Only one octant need to be considered due to the 64^ symmetry, 
so we take and Uy positive, and Ux > Uy. The segment under study starts out from q — 
(0,0). It is obvious that only coprime q — {nx,ny) can be reached. Now let VLq be the part of 
phase space (0, a) for which the trajectery hit disk q. Let us make the further approximation 
that all trajectories going from (0,0) to q have the same length Iq. We then get the following 
approximation for p(A) 



ofp(A) 




(13) 





(14) 
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Figure 1: a) The Sinai billiard with definitions of of the phase space variables (j) and a. 
b) The unfolded system with free directions (corridors) indicated. 
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Figure 2: The trace of the evolution operator (unit weight, r = 0) for the BER approx- 
imation and the exphcit periodic orbit sum. Both curves has been smoothened by a 
gaussian having width a = 0.1. The disk radius is i? = 0.1. From [15]. 



where aqlfi) is given by 

ag(0) = — j d(p d{sin a) . (15) 

This approach is not restricted to the r = case. To treat the other cases we simply average a 
power of the local instability |A(0, a)\ « Rcos{a) '-^^^^ This leads to the generalized probability 
distributions 

p.(A)«5]a,(T)J(A-/,) , (16) 

9 

with 

ag(r) = - / \K{(t),a)Yd(l) d{sina) (17) 
The zeta function may now be written 

Zr{k) = 1 - / e-''=^p^(A)dA « 1 - V a,(r)e-'*='' , (18) 



3.2 Comparison with periodic orbit calculations 

The formulation of p^(A) can be further refined. The region f2g may divided into eight parts 
with s = {g, q) where g depend on the octant from which the trajectory has arrived. This eight 
octants correspond to the eight elements of the symmetry group g G Ci^, of the billiard. The zeta 
function now reads 

Z,(fc)«l-^a,(r)e-''='= , (19) 

s 

This amounts to restricting the system to the fundamental domain and study one symmetry 
subspace Ai of the original problem, see refs. jlj, |l^. For each VLg there is an orbit s periodic in 
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Figure 3: p{A) from numerical simulation, compared with the crude approximation of 
sect 3.3. The disk radius is R = 0.05. 



the fundamental domain. We choose Is to be the length of this orbit. (It may need to go through 
an intermediate disk, but Is is well defined irrespective if s is pruned or not). Other subspaces 
may also be considered by inserting group characters, see |]l5| , ^ . We will not discuss them here 
to keep things as simple as possible. 

At this point it is interesting to check if this approximate zeta function can reproduce the 
exact trace by means of the periodic orbit sum (||), see fig 2. This plot, reaching up to / = 5.5, 
involves ~ 8000 periodic orbits. The exponential proliferation of cycles means that we cannot 
go very much higher in / by explicit periodic orbit calculations. The conclusion so far is that 
the BER approximation gives a very good description, which is generally improved for larger I, 
see ref. |l^. The fine scale structure can of course not be reproduced due to the introduction 
of deltafunctions in section 3.1 and due to the BER approximation itself. The topological zeta 
function (r 1) was examined in ref. |l5|, also with good results. It is also demonstrated in ref 



15 that the exact topological zeta function is obtained in the limit R ~* 0. The conclusion is 



that the method works for any value of r. 



3.3 A crude approximations to Pr(^) 

We will now work out a very crude approximation to Pt(A) which will turn out to be both 
illustrative and useful. For any disk radius < i? < 1/2 there is a finite number of directions 
along which a trajectory may go without ever bouncing off a disk (see fig lb)! We call them free 
directions, or corridors. Consider the direction vector {x,y). A direction cannot be free if x/y is 
irrational. So we may take {x,y) as positive coprime integers {nx,ny) such that n^jriy < 1. It is 
a little confusing that {rixjUy) may refer both to a particular disk in the unfolded system and a 
corridor but thus ambiguity will turn out to be rather convenient. The direction {nx,ny) is free 
if @0 

2R < ^ . (20) 



This means that if disk q — {n,j;,ny) lies inside radius 1/2R there are no obstructing disks in 



front of it, see fig lb. Next we assume that the disk radius is small and set Iq — + n^. The 
integral (|l^ for Iq < 1/2R is trivial 

■ (21) 



q 



In order to find an approximate expression for Pr(A) we must know the density of coprime lattice 
points with respect to the radius in the first octant dc{r). This is, to leading order, found to be 
dc{r) « ^r, see Appendix and fig 6. This yields 



If? rf — )2 

pAA) ^l^SiA- lq)aq{r) = 4(A) • a,(r)|z,.A = Y^fT^^'"'^' ■ (22) 
1 ^ ' 



Obviously the disk radius has to be small for this to apply. 

If A > 1/2R the disks start the obscure each other. The accesable disks are those along the 
corridors, se fig lb. They may be written on the form g„ = q' + n-q". Here, q" = (n", n'y) is a free 
direction and q' is one its neighbours in the Farey sequence of order n" jl8|. The Og's decreases 
as H 

1 

4(^) - ^^^^ (23) 
and the corresponding lengths grow linearly 

Iq^ ~ const + n ■ qfree (24) 

The last two equations impliy that 

MA)-^ A»i^ (25) 

Of course, there is a transition region around A ^ 1/2R but if we accept a kink at A = 1/2R 
and require Pt(A) to be continous we obtain 

{4 r(^-^)^ _R-^~"'"A"'" A < 1/2R 

!r(fc^)^./.-3^-./.-. " ' (26) 
3 r(2-r) A > l/2i? 

We have multiplied the prefactors with a factor 13/12 so that Z^"ff^(0) = 0, our approximations 
otherwise violate the normalization of p(A). In fig. 3 we compare this expression with a numerical 
simulation of the billiard for the case r = 0. Our crude approximation is a good smoothened 
approximation to p(A). 

The most obvious length scale in the Sinai billiard is the length of the shortest periodic 
orbit Iq « 1. We have now learned that there is another important length scale, Iscaie = l/2i?, 
given by the kink above, being the mean free path in the unfolded system. We will subsequently 
understand that this is a very important length scale even for the quantum problem. 

A crude estimation of the zcta function Z^'""'^'^(fc) is obtained by a Fourier transform of 
pcrude^^y What Can we say about its analytic structure? 

The leading zero fco = —ihij) as obtained from this zeta function is plotted versus r for two 
different disk radii in fig. 4. We see that the behaviour is not at all linear as was the case in the 
example of section 2.2, which is the case only for uniformly hyperbolic systems. The Sinai billiard 
is anything but unifomly hyperbolic, the intermittency becomes more and more pronounced as 
R decreases and the deviation from a linear function increases accordingly. Note that a linear 
dependence is often assumed in some oversimplified (erroneous) estimations of the entropy harrier 
in the litterature. 

For small T,that is, when the zero is close to the origin, the zeta function is dominated by the 
large I tails. Further away from the origin it is dominated by the closest disk. The topological 
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Figure 4: The generalized entropies /i(t) as given by the leading zero of the crude ap- 
proximation of the zeta function. 

entropy h{l) obtained this way is 5% wrong in the limit R ^ 0, see ref. The reason is 

twofold: the ad hoc renormalization of the prefactor, and secondly, for this large r the value is 
dominated by the closest disks. 

Equipped with our crude estimation of the zeta function we can insert it onto (|^) and the 
compare the trace obtained with the result of section 3.2. The comparison is made in fig. 5. Wee 
se that the crude approximation give a good average description. We see that when I > Iscaie the 
curve starts to approach an asymptotic behaviour. What is this limiting behaviour. 

The power law tail of Pr(A) introduces a branch cut of the zeta function along the positive 
imaginary axis. The asymptotics of the Fourier transform (|^) is governed by this cut. Let us be 
slightly more general and consider an arbitrary (generalized) distribution Pr (A) with a power law 
tail ~ 1/A™. Then an elementary calculation yields 



Aq is any point in the tail of Pr(A). The results are valid also for non integer ni. To obtain these 
results we let the contour of the Fourier integral run along the real axis. To get consistency with 
the explicit periodic orbit sums (|) for T > one has to extend the integration below all the zeros 
which will pick up up exponentially increasing terms. 

Note that there will be a sudden change in the asymptotic behaviour when t = 0. One say 
that a phase transition occurs at this point; the exponential behaviour of the trace is taken over 
by a power law. When t — Q the cut reaches down to the leading zero preventing exponential 
mixing. 

3.4 Relation to cycle expansions 

All values of the symbol s — {q,g) (introduced in section 3.2 and corresponding to one laminar 
segment) can be realized if q lies within the horizon. An idea is now to build strings of these 




Z{0) = m > 2 
Z{0) 7^ m > 1 



(27) 
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Figure 5: The trace of the evohition operator (r = 0), calculated from the crude approx- 
imation (sect. 3.3) and the more refined approximations (sect. 3.2) of the zeta functions, 
for disk radius i? = 0.1. 
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Figure 6: The number of coprime lattice point, with and without Maslov weigths, inside 
radius r, compared with the asymptotic results obtained in the Appendix. 
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symbols and build up arbitrary trajectories S1S2S3 ■ • ■ and so define a symbolic dynamics. If all the 
s's lies inside the horizon almost all such sequences could be realized and the symbolic dynamics 
would be almost complete. This would yield a fundamental part of a cycle expansion 

If we compare this expression with eq (18) together with cq ( ^l| ) we see an apparent similarity. 
Indeed the prefactors are almost the same, see (To make this comparison one first writes 

the sum J2s={q.,g) = Eg Eg and sums over g.) 

The cycle expansion breaks down when some of the s's are outside the horizon whereas the 
BER approximation works fine. The conclusion is that the approximation worked out in sections 
3.1 to 3.3 provides a good generahzation for the fundamental part of a cycle expansion. Indeed it 
is both better (it automatically preserves unitarity) and simpler as it does not use an enourmous 
amount of periodic orbits to explore the available phase space for one single symbol but rather 
measure it directly. Neither does does it need all these cycles in order to find out that the invariant 
density is uniform as we already knew that in advance. 



4 The Sinai billiard - Semiclassical considerations 

For the quantum Sinai billiard a major part is played by the (one parameter families) of neutral 
periodic orbits. The level density can be written as a sum of the mean density and two oscillating 
parts, one from the unstable orbits (expressed by the Gutzwiller formula) and one due to neutral 
orbits. 

diE) ^Y.^iE- E^) ^ d{E) + dZ{E) + d°riE). (29) 

V 

The semiclassical approximation can be refined by taking creeping orbits jl^ and ghost orbits 
into account , but we neglect them here. We we will now discuss the two oscillating contributions 
separately. In the discussion of the contribution from unstable orbits we will focuse on convergence 
and singularities of the zeta function. The considerations for neutral orbits are worked out for 
use in section 5. 



4.1 The contribution d°^'^{E) from unstable orbits 

We use the BER approximation to approximate the Gutzwiller- Voros zeta function and we will 
use the rough approximations of section 3.3, putting t — 1/2, and k = V2E = —k. According 
to our previous investigations we now expect this zeta function to approximate the Gutzwiller- 
Voros zeta function for small k, well below the first quantum state. This is certainly not the 
semiclassical limit but will lead to some interesting observations. 

We now have to introduce Maslov indices. The generalization oip^- now includes an oscillating 
weight (we still restrict us to the Ai subspace): 

Pi/2(A)« E (-l)"^+"''+'a,(r)5(A-?g) . (30) 

The aq(T)'s are the same as before but the definition of the generalized density of coprimes dc{r) 
must now include the oscillating weight. The result appears to be exactly one third of the previous 
result dc ~ If '"^i Appendix and fig 6. The crude approximation of the probability distribution 
is obtained by simply multiplying eq. (|2^) by one third. 

The zeta function is obtained by Fourier transforming our approximation of pi/2(A). An 
interesting question is if its leading zero is still below the real k axis. In this approximation it 
follows from (^) (with the extra factor of one third) that this is the case if 

i?<^i;i^. 0.3905 (31) 
135y2r(3/2) 
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For radius as big as the limiting value the approximation has already ceased to be relevant. We 
must expect that the Euler product @) is divergent on the real axis and in a strip below, but the 
leading zero can be very close to the real axis. It may very well happen that the zero crosses the 
real axis for other subspaces. 

If there is a zero below the real axis it will mean that the Gutzwiller trace formula diverges 
on the real axis. In that case, whenever we encounter as periodic orbit sum like it must be 
regularized. This amounts to keep the contour of the Fourier integral along the real axis. As 
all our approximate zeta functions converges in the entire lower halfplane this is not a major 
problem. 

It is interesting to compare this situation with the Selberg zeta function for systems with 
constant negative curvature. Along the real fc(= —k) axis the Selberg zeta function has zeros 
on the exact (quantum) locations, which is a peculiarity of these systems. In addition it has a 
leading zero on the negative imaginary axis, in exact analogy with our crude estimations of the 
zeta functions for the Sinai billiard. The situation is similar for the Riemann zeta function. 

Our semiclassical zeta function has a cut along the positive imaginary axis (and an infinite 
number of parallell cuts, cf. ref. Do we allow the semiclassical zeta function to have non 

analytical features? Well, the semiclassical approximation to the spectral determinant D{k) is 
D{k) « Dsc{k) = exp(— i7r(iV + N°'"'))Zgv{k)- Its exact quantum counterpart D{k) obeys the 
functional equation D(k) = D{—k). There is nothing saying a priori that Dsc{k) must obey this 
relation (it would be violated by a cut), it could be an artifact reflecting the shortcomings of the 
stationary phase approximation and/or the singularity of the semiclassical limit. 

But it should be noted that the best semiclassical calculations for chaotic systems so far is 
obtained by imposing the functional equation to the Gutzwiller- Voros zeta function. It is shown 
I21I that this requires a remarkable bootstrapping property between short and long periodic 
orbits. That such a bootstrapping indeed occurs is supported by the results of ref ||22| . There 
are ways of saving the semiclassical determinant from violating the functional equation. The 
contribution from neutral orbits has cuts at the same places as Zqv but they do not cancel, but 
the mean staircase N{k) could have nonanalyticities cancelling the (main) cut of Zqv and N°'"^ . 
This would extend the bootstrapping ideas mentioned above as N is governed by the short time 
behaviour and the cuts of Zqv by the long periodic orbits. Generally N is given by an asymptotic 
expansion so it shouldn't be surprising if this expansion lead to a cut. 

In the classical case the cut carry information about the power law decay of correlations [|l5j . 
It would be nice to know what it tells us in the quantum case. 

4.2 The contribution d'^'^{E) from neutral orbits 

There is a neutral orbit in any free direction. The general formula is [^3| 

+ 

dl'" = '^2BiCos{Kk) (32) 

with 



(27r)-^/^ y \ni\K. 

The index i runs over all neutral periodic orbits. Their length is written as li = riilp. where is 
a repetition number and /p. the length of the corresponding primitive periodic orbit. Dp- is the 
geometrical width of the orbit. Note that this contribution is 0{1/E^^'^) whereas the one due to 
unstable orbtits is Oil/E^/"^). The reason that the unstable orbits still are important is because 
they are so numerous. Generally the neutral orbits affect the low part of the spectrum and large 
scale structures in the higher part of the spectrum. 

For later puposes it is convenient to sum over both positive and negative traversals 



gi«".ip. (34) 
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Instead of applying these equations directly it is instructive to first consider the integrable R ^ 
limit of the Sinai billiard. The wavefunction must now vanish at the midpoint which is the relic 
of the disk. The spectrum is given by _E = 27r^(n^ + m^) where m,n are postive integers. This 
leads to a mean level density d = I/Stt which differs form the Weyl term d = 1/2-k. The reason is 
that is that the i? is a very odd system quantum mechanically, only having wave functions 
in two of the symmetry classes {A2 and i?2). 

To find d'^"'^{E) we perform the Berry- Tabor trick and rewrite the level density by means 
of the Poisson's summation formula 



oc 

d^{E) =^S{E~ 27r2(TO2 + n^y^ ^ ^ I dn dm 6{E - 2TT^{m^ + „2))g2x^(A^„+Mr„) (35^ 

m.n M,N=-oo 

Performing the integral by stationary phase we arrive at (^2|) with 

1 

(27r)3/2V7V2 _^ M2 

and 



B^ = ^ (36) 



h = VW+M^ (37) 

where i ranges over all pairs {AI,N) in the first quadrant. Only one symmetry class (say A2) 
is considered if we restrict the summation to the first octant (not bothering how to treat lattice 
points on the symmetry line M — N). If we now write the lattice vector as (M, N) — n{nx, Uy) 

where {nx,ny) are coprime, we will get consistency with ( p3| ) if Ip — \ /n'^ + and Dp = l/lp. 

The width Dp decreases as the radius R increases. Inspectin g (pO[ ) we realize that the width 
then becomes Dp = l/lp~ 2R. These results are derived in ref. [[17| in another way. 

With the equipment procured so far we will now attempt to calculate the spectral form factor. 

5 More semiclassics - The spectral form factor 

A lot of information about spectral fluctuations is encoded in the spectral form factor 

K{T) = ^ / dee''^'^'^'{d°'^{E - -)d°''{E + -))e ■ (38) 

2 2 

The averaging {■■■)£ has been performed in a small range AE such that d^^ ^ AE ^ E. If 
there are contributions from both neutral and unstable orbits the form factor divides into three 
parts (with obvious notations) 

K{T) = if„„(r) + KU^) + i^„„(r) (39) 

The suggested universal form factor relevant for time reversable chaotic systems is that of the 
Gaussian Orthogonal Ensamble (GOE) which reads ||2^ 

r 2r-Tiog(i + 2r) r<i 
^«o^ = \ 2-Tiog|f±i r>i 

5.1 The contribution from neutral orbits: Knn{T) 



If we insert d'^^ into the definition of iir„„(T), and expanding ^/2E + e « v2i? + t/^2E we 
obtain ^ 

^-(^) = f EE^»^.'5(^-^-^)w^a^-^.-))>i^ • (41) 

i 3 



14 



Unstable+Neut ral 
Unstable 
Neutral 



0.2 

tau 



Figure 7: The contributions Knn and K^u and their sum. The latter is calculated from 
eq (44) using the methods of section 3.3. The parameters are R = 0.1 and Tgcaie = 0.1. 



where the length variable Z(T) = ^f^E^ndT . This expression is always dominated by the diagonal 
terms il 

K^^{r)^^Y.^l8{l-U) . (42) 

i 

For an integrable system this always equals unity [|4|. This can be easily checked for the Sinai 
R billiard by replacing the sum J2m n integral and using polar coordinates. 

If R is very small and / is far from Iscaie-, Knn will still be close to unity, but it will now 
decrease with increasing T. An example is given in fig. 7. 

Let us now study the case I ^ I scale- We can be more general and consider and arbitrary 
system with a finite number of primitive neutral orbits in the limit / 3> maxp lp{ = I scale for the 
Sinai case) we have 

- (^-E^.^IE^('-w = (^-^E^.^^ ■ (43) 

We thus have an inverse power law decay, Knn ~ contradiction to the large T limit of 

the GOE prediction Kqo e ^ ^ — iW^ ■ 

5.2 The contribution from unstable orbits: KuuC^) 

In the limit T <^ 1 this contribution to the form factor is also dominated by the diagonal terms. 
However, the result must be multiplied by two |25[ since most periodic orbits in time reversable 
systems are rotations, and each rotation counts as two cycles with exactly the same length. We 
insert the trace formula into the definition of the form factor and express the result in terms of 
the trace of the evolution operator 

Knu{r) = 2Ttr£'}^ (44) 
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Figure 8: The spectral rigidity of the Sinai bihiard calculated from two subsets of the 
spectrum and compared with GOE. The radius R = 0.1. The relevant scale Tgcaie is 
slightly above (below) Tg^aie = 0.1 for the lower (higher) part of the spectrum. 



An example of this function is found in fig 7. If Tscaie T <C 1 we find, after consulting eqs. (Em 
and ( p7[ ) , that K^u = 2T — 2Tscaie where the first term agrees with the universal GOE result p5[ 
and the second being a correction dying out in the limit iJ — > oo. Tscaie is the scale corresponding 

to Iscale given by %cale = lscale/{V2E2TTd). 

The big challenge for semiclassical methods is to calculate the large T limit of Kuu{T). 
General considerations j25j say that the form factor must approach unity as T ^ cx3. The 
statement is often written in terms of an explicit periodic orbit double sum analogous to (p[). 
One should always be extremely careful when encountering such sums since they in general require 
regularization even if the double sum is convergent, otherways one will obtain an exponentially 
increasing form factor |^ . 



We will not calculate the mixed contribution Kun', it may provide small corrections. It seems 
unlikely that K^n and/or Kuu could cancel the 1/T decay of Knn- Summing up our contributions, 
fig 7, we conclude that for any finite energy the form factor is not compatible with GOE for small 
T either. The small T limit governs the large scale structure of the spectrum. A convenient 
mesaure to study such properties is the spectral rigidity which we now will discuss. 



5.3 Spectral rigidity 

Spectral rigidity, or A3(L) statistics measures the (mean square) deviation of the spectral 
staircase from a linear function over L mean spacings. Small values of A^i^L) thus correspond to 
a rigid spectrum. A3(L) is expressed in terms of the form factor as p5| 

^'^^^ =2^1 (45) 
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where 

^, , 1 ^cos'^(y) ^cos(y)sin(y) „sin^(y) 

y2 yZ yi yi 

The function G{y) is close to unity if y > 1 so for sufficiently big L the rigidity explores the small 
T limit of the form factor. 

The deviation from GOE for small T implies that even if A^{L) follows GOE for a while it 
must certainly deviate from it well before L — l/7^caZe = '2RV'2E2T:d. To study this we have 
also computed the exact quantum spectrum by the methods of ref. |l7|. It is clearly seen in fig. 
8 that such deviation occurs well before the saturation value Lq corresponding to the shortest 
periodic orbit Iq- Such behavior has been observed before |^. If the form factor approaches a 
constant c as T — > 0, then the rigidity will be asymptotically linearly increasing with slope c/15. 
The slope in fig. 8 is roughly 0.17 corresponding to a constant c ^ 0.25, which is quite compatible 
with the results of fig 7. As we said, we expect c to approach unity when the disk radius R 0. 
Accordingly the slope of A3 should approach l/15_(as for an integrable system). That this is 
indeed the case is supported by the results of ref. Q . 

The effect we have seen is mostly due to the neutral orbits. However, we can expect similar 
deviations, but perhaps less pronounced, for systems without neutral orbits. The reason is that 
the asymptotics of the trace (which was a slow power law for the Sinai billiard) does not set 
in until a length scale Iscaie which can be much longer than the shortest periodic orbit Iq- For 
instance, in the hyperbola bilhard Iq = I whereas Iscaie ~ 27.61, see ref. [||. 



6 Concluding remarks 

We have surely not used the full potential of the methods outlined in section 3. I conclude by 
listing some lines along which the work may proceed. 

• Improve the BER approximation by replacing the delta functions in sec 3.1 by something 
more realistic. 

• Find corrections (similar to curvature terms) to the BER approximation. One should also 
correct for the nonmultiplicativity of the weight. 

• Try to calculate the long T limit of the form factor. It is thus essential to have good control 
of the fine scale structure the (properly regularized) trace trCl^^^^^- requires that the 
previous two points have been worked out. Is there perhaps any simple universal structure 
in this fine structure? 

• Work out the BER approximation (with corrections) for other systems. Each system will 
of course present its own problems: the stadium billiard has no natural division between a 
chaotic and a regular part, as it is glued together from regular components. The anistropic 
Kepler problem and the closed three disk billiard suffers from infinite sequences of 
cycles accumulating to finite length and stability. The hyperbola billiard has infinite horns. 

This work was supported by the Swedish Natural Science Research Council (NFR) under 
contract no. F-FU 06420-303. I will also thank Predrag Cvitanovic for various remarks. 

Appendix 

The number of coprime lattice points inside radius R in the first octant, Nc{R) 



N,ir)= J2 Q{r-^nl+nl) (47) 

Scd(nj;,ny) = l 



is given by the integral equation 

_ NJm = / - 

dr r 



(48) 
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This equation expresses the fact that density of coprimes at radius r contributes to the density 
of non-coprimes ar radii 2r, 3r etc. Each non-coprime is uniquely related to one coprime point. 
This summability leads to a simple integral as above. The total number of lattice points inside 
radius R is ttB? /S, appearing in the left hand side, [x] denotes the integral part of x, which is on 
the average = x —1/2. Inserting this into the equation above it is easily seen that Nc{r) = j^r'^ 
is a solution. 

The introduction of maslov indices in section 4.1 motivates us to study the asymptotic be- 
haviour of 



Numerical calculation suggests that iVc(r) ~ 1/3 • Nc{r). How can this be understood? The pairs 
{nxiUy) can either be (odd,odd), (odd, even), (even, odd). The (even, even) option is excluded by 
not being coprime. If the three possibilities have the same density the result follows. 
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